Lattice Vibrational Modes in Si/Ge Core-shell Nanowires 
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We present a first-principles study on lattice vibrational modes of Si/Ge core-shell nanowires 
(NWs). In addition to quantum confinement, the internal strain induced by the lattice mismatch 
between core and shell contributes to significant frequency shifts of characteristic optical modes. 
More importantly, our simulation shows that these frequency shifts can be detected by Raman 
scattering experiments, providing convenient and nondestructive ways to obtain structural informa- 
tion of core-shell materials. Meanwhile, another type of collective modes, radial breathing modes 
(RBMs), are identified in Si-core/Ge-shell NWs and their frequency dependence is explained by an 
elastic media model. 



I. INTRODUCTION 

Core-shell nanowires (NWs), a type of novel quasi-one- 
dimensional nanostructures, have ignited tremendous in- 
terest to date because these unique radial heterojunc- 
tions own extra degrees of freedom to tune the mechan- 
ical, electric and optical properties by varying their core 
and shell sizes, respectively, which make them superior to 
usual homogenous materials^—. Among numerous core- 
shell structures, Si/Ge core-shell NWs are of particular 
interest because of the known importance of silicon on 
microelectronics. Recent research has shown that Si/Ge 
core-shell NWs own a number of unusual thermal trans- 
port features for energy applications^—. For example, 
the enhanced difference in thermal and electric trans- 
port behaviors between the core and shell gives rise to 
higher-efficiency thermoelectric devices based on core- 
shell NWs£. 

In order to realize and optimize the above interesting 
properties, it is imperative to determine how to obtain 
structural information and, thereafter, how to design core 
and shell geometries. Traditional approaches, such as 
the transmission electron microscopy (TEM) and neutron 
scattering measurements, are time consuming or expen- 
sive for the large-scale production; alternative convenient 
ways are highly desired. On the other hand, lattice vibra- 
tional modes are known to be sensitive to the chemical 
bonding and boundary conditions of nanostructures and 
they are well accessible by varieties of experiments, such 
as Raman scatterin g 9 ! 10 ' 12 , making it possible to study 
lattice vibrational modes of core-shell NWs and extract 
useful structural information. 

Beyond application interests, lattice vibrational modes 
and corresponding phonons are fundamental properties 
of solids. However, compared with bulk semiconduc- 
tors and homogenous nanostructures, there have been 
very limited first-principles studies on vibrational modes 
of core-shell NWs although the core-shell structure will 
surely introduce richer physics. Therefore, calculating 
lattice vibrational modes in Si/Ge core-shell NWs will be 
an appropriate starting point to understand electric and 
thermal properties of reduced dimensional structures. 

In this work, we employ first-principles calculations 
and focus on the following two topics of Si/Ge core-shell 



NWs: 1) how the unique core-shell structure modifies lat- 
tice vibrational modes and what their impacts on ther- 
mal properties are; 2) how to extract the structural infor- 
mation from characters of lattice vibrational modes and 
provide clues for experimental measurements. 

Through our simulations, we find that the structural 
information of core-shell NWs can be efficiently obtained 
by frequency shifts of high-energy optical vibrational 
modes. This type of frequency shifts can be well ex- 
plained by the change of the strain induced by the lat- 
tice mismatch between the core and shell. Moreover, we 
have calculated the Raman scattering spectrum of Si/Ge 
core-shell NWs, showing that these frequency shifts can 
be detected by experiments. Meanwhile, we identify the 
existence of RBMs in Si-core/Ge-shell NWs, whose fre- 
quency variation can be understood by an elastic media 
model. 

The rest of this article is organized as the following. 
In Sec. mi we present our calculation methodology and 
computational setup. In Sec. IIII1 we discuss the simula- 
tion results, including properties of high-frequency opti- 
cal modes, Raman scattering spectra, and the RBM fre- 
quency dependence on the structure of core-shell NWs. 
In Sec. IIV1 we summarize our studies and conclusion. 



II. COMPUTATIONAL DETAILS 

Our calculation is based on the density functional 
theory (DFT) within the local density approximation 
(LDA)ii^ as implemented in the Quantum ESPRESSO 
package^. We employ the plan-wave basis and the pseu- 
dopotential approximation^. The plane-wave energy 
cutoff is 16 Ry. For the Brillouin zone (BZ) integra- 
tion, we use a 1x1x8 k-point sampling grid. All cal- 
culations are done in a supercell arrangement^ with a 
0.7-nm vacuum distance between neighboring NWs. Our 
studied NWs are oriented along the [110] direction, an 
energetically preferred configuration^ 9 -. Those dangling 
bonds on the surface of NWs are passivated by hydrogen 
atoms. All NWs are fully relaxed according to the force 
and stress calculated under DFT/LDA. 

There are a number of degrees of freedom to decide the 
structure of core-shell NWs, such as the core size, shell 
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TABLE I: The structure information of Si/Ge core-shell NWs and homogenous NWs investigated. (SiGe means Si core and Ge 
shell while GeSi means Ge core and Si shell.) 



NWs 


Core 


Shell 


Core diameter (nm) 


Shell diameter (nm) 


Lattice constant along wire (nm) 


ol-l 


Ol 


ol 


9 9/1 


9 9/1 


U.o / o 


GeSi-l 


Ge 


Si 


1.12 


2.21 


0.383 


GeSi-2 


Ge 


Si 


1.75 


2.30 


0.390 


Ge-1 


Ge 


Ge 


2.41 


2.41 


0.401 


SiGe-1 


Si 


Ge 


1.09 


2.33 


0.396 


SiGe-2 


Si 


Ge 


1.66 


2.29 


0.388 


Si-2 


Si 


Si 


1.61 


1.61 


0.382 


Ge-2 


Ge 


Ge 


1.76 


1.76 


0.405 



size and the chemical compositions. Because it is impos- 
sible to study all of them in one article, we investigate 
Si/Ge core-shell NWs with a fixed diameter (~ 2.3 nm) 
but a varying core-shell size ratio. At the same time, ho- 
mogenous NWs are also studied for comparison purposes. 
The detailed structural information of our studied NWs 
is reported in Table Q] 

We employ the linear response approach to obtain lat- 
tice vibrational modes and their frequencies 2 ^ - — . Only 
vibrational modes at the T point of the first BZ is calcu- 
lated because it is of particular interest for the first-order 
Raman scattering measurement. 

We calculate the Raman scattering spectrum by con- 
sidering the non-resonant first-order process, which can 
be described within the Placzek approximatio n ' 25 . 
Then the Raman activity l k amman associated with the 
vibrational mode k is 

Iramman l e s ' Oi • (1) 



where e s is the polarization of incoming photon and ex, 
is the polarization of scattered outgoing photon. The 
Raman tensor a k is defined as 



" fe = ^Ef7^, (2) 

, ul al 

at 

where \ij is the electric polarizability tensor, a is the 
index to specify the atom in the unit cell, and I is the 
index of coordinates, u^, is the displacement of the atom 
a along I direction in the vibrational eigenmode k. 

For realistic experimental cases, NWs are usually ori- 
ented randomly. Because of the depolarization effect, we 
may only consider the circumstances that the direction 
of the incident beam, the polarization direction of this 
beam, and the direction of the observation are perpendic- 
ular to each other. Then the spatially averaged Raman 
activity of those randomly orientated NWs is given by: 



Iramman = 5(«4 + <4v + a zzf + j[("L - ^yyf + («L ~ <4 t ) 2 + (a k yy - u\ z f + 6(\a k y \ 2 + \a h xz \ 2 + \a k yz \ 2 )}. (3) 



Because our calculated Raman spectra only consider 
the first-order process, they may be different from some 
Raman experiments which could be dominated by reso- 
nant effects. However, our obtained Raman activity, e.g., 
active or not, shall still be useful to understand those res- 
onant experiments. 



III. RESULTS AND DISCUSSIONS 
A. High-frequency optical modes 

The density of vibrational modes (DVM) at the T point 
is presented in Fig. [TJ It has to be pointed out that 
those vibrational modes mainly consisting of hydrogen- 
atom motions are excluded from our plots because we are 



interested in intrinsic vibrational modes of NWs that are 
not sensitive to the environment. Therefore, our results 
and conclusion will be more universal and less influenced 
by complicated passivations as in various experiments. 

The most prominent feature of these DVMs is that 
there are two characteristic peaks among most core-shell 
NWs, which are marked by Osi and Oc e in Fig. [TJ re- 
spectively. To understand the physical origin of these 
peaks, we have plotted the real-space motion of two typ- 
ical modes corresponding to these peaks in Fig. [51 

The first character of these modes is that they are 
mainly confined within the core or shell, respectively. In 
another word, these Si and Ge vibrational modes are de- 
coupled. For example, Fig. [2] (a) is mainly the Si-Si 
vibrational mode confined within the core regime while 
Fig. [5] (b) is mainly the Ge-Ge mode within the shell 
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FIG. 1: DVM at the T point of Si/Ge core-shell NWs and ho- 
mogenous NWs with wire diameter of about 2.3 nm in Table 
U (a) Si-1 NW (b) GeSi-1 NW (c) GeSi-2 NW (d) Ge-1 NW 
(c) SiGe-1 NW (f) SiGe-2 NW. The red dashed line is used 
to guide the shift trend of the peak of high-frequency Ge-Ge 
optical modes. 



regime. This decoupling is due to the significant mass 
difference between Si and Ge atoms and their bonding 
strength differences. 

The second character is that the neighboring atoms 
of these modes plotted in Fig. [2] have opposite vibra- 
tional directions, a typical feature of optical modes of 
semiconductors at the T point. Based on the above two 
characters, we conclude that these peaks (Osi and Oce) 
originate from those high-frequency optical modes of bulk 
Si and Ge. This can be further confirmed by their fre- 
quency; the Oce peak has a frequency around 8 THz 
while the Osi peak has a frequency around 15 THz. They 
are consistent with frequencies of optical modes in Si and 
Ge bulk crystals. 

However, the frequencies of peaks (Oce and Ost) are 
not exactly the same as their bulk counterparts. Instead 
they exhibit substantial frequency shifts as the core-shell 
size varies as shown in Fig. [1] For example, as we expand 
the Ge-core size from zero to the whole Ge-core/Si-shell 
wire (Figs. [T] (a) to (d)), both peaks Osi and Oc e have 
a red shift. Similarly, both peaks have a blue shift as Si 
core size increases in the case of Si-core/Ge-shell NWs, 
as shown in Figs. d](d) to (f). 

The physical origin of the above shifts of prominent 




FIG. 2: The typical vibrations of highest-frequency optical 
modes in Si/Ge core-shell NWs. Hydrogen atoms on the sur- 
face are not plotted. The arrows stand for the motion di- 
rection and magnitude, (a) The highest-frequency Si optical 
mode in the Osi peak of a Si-core/Ge-shell NW. (b) The 
highest-frequency Ge optical mode in the Oae peak of the 
Si-core/Ge-shell NW. 



peaks in the DVM can be attributed to two factors, quan- 
tum confinement and strain condition. First, quantum 
confinement in such narrow NWs definitely modifies the 
frequency of lattice vibrational modes as discovered by 
previous studies^—. However, since the whole diameter 
of our studied NWs is fixed, quantum confinement shall 
not be a major factor for the frequency shift observed in 

Fig. m 

On the other hand, the internal strain of core-shell 
NWs induced by lattice mismatch can be an important 
reason for those frequency shifts. In our Si/Ge core-shell 
NWs, the Si part is stretched while the Ge part is com- 
pressed. As the core-shell size changes, the average lat- 
tice constant along the axial direction varies according to 
the Vegard's law (a = a° A {\ — X) + a B X , where a is the 
lattice parameter of A\_x Bx crystal, X is the concentra- 
tion of constituent element B, and a B are the lattice 
parameters of pure A and pure B crystals, respectively), 
which can be seen from Table I. As a result, the intrin- 
sic strain will gradually change, resulting in significant 
frequency shifts of Si and Ge modes. 

To further confirm the above explanation, we focus 
on the frequency shifts of the particular optical modes 
picked from Osi and Oce peaks, respectively. As shown 
in Fig. |31 for both Si-core/Ge-shell and Ge-core/Si-shell 
NWs, these frequency shifts are approximately linear to 
the square of the ratio of the core and NW radii, which 
is also proportional to the ratio of the number of atoms 
in the core and the whole NW. Since this linear relation 
is exactly the lattice constant variation according to the 
Vegard's Law, the change of the frequency shall be a 
result of the variation of the internal strain. 

Meanwhile, we calculate homogenous Si NWs (SiNWs) 
and Ge NWs (GeNWs) with applied axial strain to check 
the strain effect as shown in Fig. HJ The size of homoge- 
nous NWs is chosen to be similar to the corresponding 
core size of Si/Ge core-shell NWs (SiGe-2 and GeSi-2 
NWs in Table I) and we focus on the highest-frequency 
optical mode, which is the same mode plotted in Fig. 
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FIG. 3: The highest frequency of the optical mode picked 
from Osi and Oce peaks, r is the core radius and R is the 
whole wire radius, (a) Si-core/Ge-shell NWs with a fixed 
wire diameter, (b) Ge-core/Si-shell NWs with a fixed wire 
diameter. 



FIG. 4: Highest optical mode frequencies with external uni- 
axial strains applied to pure Si NW and pure Ge NW respec- 
tively, (a) Si-2 NW, the cross marks the frequency position of 
the highest Si optical mode in SiGe-2 NW. (b) Ge-2 NW, the 
cross marks the frequency position of the highest Ge optical 
mode in GeSi-2 NW. 



[3 For both SiNW and GeNWs, these modes exhibit 
a monotonic change with the applied axial strain; their 
frequencies increase with the compressive strain and de- 
crease with the tensile strain, which is consistent with the 
trend of the peak shifts in core-shell NWs. Moreover, we 
mark the highest-frequency optical mode of the corre- 
sponding core-shell NWs in Fig. 01 We find that their 
frequencies are close to those of the highest-frequency op- 
tical mode of homogenous NWs under the similar strain. 
Therefore, the variation of the strain condition is the 
main reason for the frequency shifts of high-frequency 
optical modes in core-shell NWs. This result builds a 
bridge to connect the frequency shift of characteristic op- 
tical modes with the structure of core-shell NWs. 



B. Raman scattering spectrum 

Raman scattering is a widely used approach to detect 
lattice vibrational modes of solids. Recently it is calcu- 
lated in SiNWs as wel l 12 ' 13 . If the above frequency shifts 
of characteristic optical modes can be observed in their 
Raman scattering spectrum, it will be useful to under- 
stand the structures of core-shell NWs. Following the for- 
mulas in Eq. ([2]) and Eq. ([3]) , we calculate the first-order 
Raman scattering spectra of our studied Si/Ge core-shell 
NWs and present them in Fig. [SJ 

In Fig. [SJ we can see the following important fea- 
tures. First, those highest-frequency optical modes in 
the peaks Osi and Oce are strongly active in the Ra- 
man scattering spectra. This is reasonable because they 
originate from bulk optical modes that are highly Raman 
active. Since we have shown that these shifts are caused 
by the strain and associated with structural variations, 
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FIG. 5: Raman activity of Si/Ge core-shell NWs and homoge- 
nous NWs. (a) SiGe-1 NW. (b) SiGe-2 NW. (c) Si-1 NW. (d) 
GeSi-1 NW. (e) GeSi-2 NW. (f) Ge-1 NW. Osi marks the 
highest-frequency Si optical mode, Oge marks the highest- 
frequency Ge optical mode, and B marks the RBM. 



the active Raman scattering peaks provide a convenient 
way to obtain structural information and corresponding 
strain condition inside Si/Ge core-shell NWs. As a re- 
sult, they will be able to answer an fundamental question 
in nanoscience, which is if the huge strain can sustain in 
narrow core-shell NWs without introducing an significant 
amount of dislocations. 

We also observe that the Raman scattering signal in- 
tensity strongly depends on the type of atoms involved 
in the vibrational modes. As shown in Fig. [SJ those Ge 
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modes have much stronger Raman signals than those Si 
modes. This can be attributed to the larger size of Ge 
atoms, so that their motions can induce much stronger 
changes of polarizability, leading to enhanced Raman sig- 
nals. 

Meanwhile, the intensity of Raman scattering signals 
also depends on the spatial location of those atoms in- 
volved in the vibrational modes. The Raman scattering 
peak is usually enhanced if the corresponding mode is 
within the core region while it is depressed if the mode is 
within the shell region. For example, though SiGe-2 (Si 
core) NW has a lower proportion of Si atoms than GeSi- 
1 (Si shell) NW as listed in Table I, the relative Raman 
signal of the highest-frequency Si mode in SiGe-2 NW is 
stronger than that in GeSi-1 NW, as shown in Figs. [5] 
(b) and (d). The similar result can be observed in other 
NWs if we compare Figs. [5] (a) and (e). 

Other than core-shell NWs, if comparing the Raman 
spectra of the homogenous SiNW and GeNW with a 
similar diameter as shown in Figs. [5] (c) and (f), we 
find that the SiNW, different from the GeNW which 
has only one dominant Raman peak at highest opti- 
cal frequency, displays bright Raman activities for many 
phonon modes with lower frequencies. Since both bulk 
Si and Ge have only one dominant peak in the first- 
order Raman scatter spectrum, those newly active low- 
frequency modes in SiNWs imply that these vibrational 
modes are more affected by the quantum confinement 
and symmetry breaking in SiNWs. This has been ob- 
served by recent experiments^^. Such a different varia- 
tion of vibrational modes may be helpful to explain why 
the thermal conductivity of SiNWs differs from their bulk 
crystals more than that of GeNWs^. 



C. RBM 

After discussing the high-frequency optical vibrational 
modes in Si/Ge core-shell NWs, we turn to another type 
of collective modes: the RBM. Because of the unique 
geometry of quasi-one-dimensional nanostructures, their 
RBMs arc of great interests and importance. For exam- 
ple, they are characteristic modes reflecting the geome- 
try of carbon nanotube o 27 ' 28 . Recent researc h 12 ' 30 has 
also shown that the RBM are Raman active in narrow 
semiconducting NWs and they are tightly related to the 
diameter of NWs. Then an obvious question is whether 
we can observe the similar RBM in core-shell NWs. 

Our first-principles result shows that the answer for 
core-shell NWs is more complicated than usual homoge- 
nous nanostructures case. For example, the RBM are 
identified in the Si-core/Ge-shell NWs studied in this 
work as shown in Fig. [B](a). Moreover, we find that this 
RBM are Raman active and can be detected by the Ra- 
man scattering experiments, as marked in Fig. [5] (a), (b) 



and (c). On the other hand, for the other type of NWs, 
Ge-core/Si-shell NWs, we cannot identify the RBM sim- 
ply by sight. In fact, even for GeNWs with a similar 



(a) 










M |V 








FIG. 6: (a) A typical RBM in Si/Ge core-shell NWs. (b) The 
frequency of RBM according to the core/wire size ratio. The 
fitting curve is based on Eq. (JSj . 



diameter as the studied core-shell NWs, no RBM can be 
easily identified. The relative softer bonding between Ge 
atoms and the relatively larger mass of Ge atoms sub- 
stantially reduce the frequency of the RBM. As a result, 
the RBM will mix with other modes because of the ap- 
plication of the sum rule in our simulations. 

The most important motivation to study the RBM is 
to see if they have a strong relation with the structure 
of nanomaterials, such as the diameter of NWs. Pre- 
vious studies^ 2 - have successfully applied an elastic me- 
dia model to study the RBM in SiNWs, encouraging 
us to apply a similar approach to study the RBM in 
Si-core/Ge-shell NWs. By using the classic wave equa- 
tion with a cylindrical boundary condition and regard- 
ing core and shell regions as two different homogenous 
elastic medias with different sound velocities ci and Ci 
respectively, the vibration of RBM can be described as 
u(x,y,z) — u(x,y) = u(r)f and satisfies the following 
equations: 



d 2 u(r) 1 du(r) lu 



dr 2 



dr 



for < r < Ri, 
(4a) 



d 2 u(r) 1 du(r) ,uj 2 1 . . . „ „ „ „ 
' + T^ 1 + (— - — Mr) =0forR 1 <r< R 2 , 



dr 2 



r dr 



(4b) 

where w is the frequency of RBM. i?i is the core radius 
and i?2 is the shell radius, u\(r) and U2(r) are RBM 
wave function in the core and shell region, respectively. 
We solve Eqs. (|4]) using the free boundary condition at 
r = i?2 and continuous and differentiable continuous con- 
dition at r = Ri, we get the equation Eq. ([3]) which de- 
termine u), the frequency of RBM: 
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Cl 
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''2 



'-2 



'2 



(5) 



r 



where Ji and Yi are first-order Bessel functions of the 
first kind and second kind, respectively. By applying 
perturbation theory and the approximation of Bessel's 
functions in Eq. ([6]) when x | (it is satisfied in our 
simulations), 



JiO) fa J—cos(x- jn), 

V 7TX 4 



(6a) 



in Fig. EJb). From the fitting curve, the ratio of phonon 
group velocities in Si core and Ge shell is about f.82, 
which is larger than the ratio in homogenous NWs or 
bulk materials (~ 1.56). This slight enlargement of the 
phonon group velocity ratio in core-shell NWs is a re- 
sult of intrinsic strain effect and it may be of interest for 
thermoelectric applications of core-shell NWs. 



IV. CONCLUSIONS 



Yi(x) 



sm(x - 



(6b) 



Then a relation about the frequency of RBMs in core- 
shell NWs is obtained: 



i?2 , Ri Ris 
w( — + ) 

C2 Cl C 2 



constant. 



(7) 



Therefore, the frequency uj of the RBM is determined 
by the ratio of core/wire size and that of the sound ve- 
locities of core and shell. The constant in Eq. ([7]) are a 
series of values of x when J[(x) — 0, according to the 
free boundary condition. If we assume the sound veloc- 
ity does not change too much, the above Eq. ([7]) can be 
approximated to 



C2 R c 2 



(8) 



From the equation Eq. ([Jj, we can easily reduce it to 
the known RBM frequency of homogenous NWs 12 : 



I? 



in- 



constant. 



(9) 



We have applied the equation Eq. ([8]) to fit the data 
of RBMs in the studied Si-core/Ge-shell NWs, as shown 



In summary, we have studied lattice vibrational modes 
of Si/Ge core-shell NWs by the first-principles approach. 
The frequency of high-energy Si-Si and Ge-Ge optical 
modes exhibits significant shifting when the core-shell 
size varies. Our analysis shows that these shifts are re- 
lated to the variation of the strain conditions inside core- 
shell NWs. At the same time, the RBM are identified 
in the Si-core/Ge-shell NWs, whose frequency strongly 
depends on the geometry of NWs, and our elastic media 
model explains well the frequency of the RBM. More- 
over, we have performed the calculation to obtain the 
first-order Raman scattering spectra of studied NWs, in 
which the shift of high-frequency optical modes and the 
RBM can be identified, providing a convenient way to 
detect the structural information of core-shell NWs. Our 
studied vibrational modes and their frequencies could be 
useful for thermoelectric applications as well. 
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